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Abstract 

The survival problem for a diffusing particle moving among ran- 
dom traps is considered. We introduce a simple argument to derive 
the quenched asymptotics of the survival probability from the Lifshitz 
tail effect for the associated operator. In particular, the upper bound 
is proved in fairly general settings and is shown to be sharp in the case 
of the Brownian motion in the Poissonian obstacles. As an application, 
we derive the quenched asymptotics for the Brownian motion in traps 
distributed according to a random perturbation of the lattice. 
Keywords: Trapping problem; random media; survival probability; Lif- 
shitz tail 
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1 Introduction and main results 

In this article, we consider a diffusing particle moving in random traps. The 
motion of the particle is given by a simple random walk or a Brownian motion 
and it is killed at a certain rate when it stays in a trap. Such a model ap- 
pears in various models in chemical physics and also has some relations to the 
quantum physics in disordered media. We refer to the papers by Havlin and 
Ben-Avraham [S] and den Hollander and Weiss [1] for reviews on this model. 

The mathematical discription of the trapping model is given by the sub- 
Markov process with generator 

= -«:A + K., (1) 

where A is the Laplacian on L'^iW^) or t^il/") and (V^,P) a nonnegative, sta- 
tionary, and ergodic random field. Heuristically, the height of corresponds 
to the rate of killing. Let us write {{Xt}t>Q, {Px}xeM.<^ or id) for the Markov 
process generated by — kA. One of the quantity of primary interest concerning 
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this process is the survival probabihty of the particle up to a fixed time t, which 
is expressed as 



expj-^ V^{Xs)ds^ 



(2) 



From this expression, we can identify the survival probability as the Feynman- 
Kac representation of a solution of the initial value problem 

dtu{t,x) = K/\u{t,x) -V^{x)u{t,x) for (t, x) G (0, oo) x M'^, 

.(o,.)-i. 

Therefore, it is natural to expect that the long time asymptotics of the sur- 
vival probability gives some information about the spectrum of H^^ around the 
ground state energy and vice versa. This idea has been made rigorous first by 
Fukushima [6j, Nakao [13j, and Pastur [Hj (with the analysis of some concrete 
examples) in the following sense: from the annealed long time asymptotics of 
the survival probability, one can derive the decay rate of the integrated density 
of states around the ground state energy. Their arguments are based on the fact 
that the Laplace transform of the integrated density of states can be expressed 
as the annealed survival probability for the process conditioned to come back 
to the starting point at time t. Therefore, the above implication follows by an 
appropriate Tauberian theorem and, since there is the corresponding Abelian 
theorem (see e.g. Kasahara [ID]), the converse is also true. 

The aim of this article is to study a relation between the quenched asymp- 
totics of Uuiit.x) and the integrated density of states. Let us start by recalling 
the notion of the integrated density of states. To define it, we assume the 
following: 

Assumption 1. In the continuous setting, belongs to the local Kato class 
-^d.ioc- (See [5] or [12] for the definition of i^(i ioc.)@ 

This assumption is sufficient to ensure that if^ is measurable in uj as an oper- 
ator. For the notion of measurability of operators, we refer to a lecture notes 
by Kirsch [TT]. (In fact, this is slightly stronger but we need this to utilize 
a uniform bound for the semigroup e~^^^ in the proof.) Under the above as- 
sumption, the integrated density of states of H^^ is defined as follows: 

A^*(A) = Jim ^E[#{A;GN;A:,,((-i?,i?)'^) <A}], * = D or N, (4) 

where \^ ^((— i?, i?)*^) (resp. ^((— i?, -R)*^)) is the k-th. smallest eigenvalue of 
if^j in {—R-i RY with the Dirichlet (resp. Neumann) boundary condition. The 
existence of the limit in the right hand side can be proved by superadditivity 
(resp. subadditivity). 

Now we state our first result. 



2 



Theorem 1. Suppose that Assumptions\J\holds and that there exists a regularly 
varying function (p with index L > such that the integrated density of states 
associated with the operator in ([T]) admits the upper bound 



A^°(A) < exp{-0(l/A)(l + o(l))} as A ^ 0. 
Then, for any fixed x G M.'^, 

P-a.s. u^(t, x) < exp {— t/?/'((ilogt)(l + o(l))} as t — >■ oo, 
where ip is the asymptotic inverse of (p. 

The foUowing assumptions are necessary only for the lower bound. 
Assumption 2. {Moment condition) There exists a > such that 



E 



sup exp{K;(x)"} 

,a;G[0,l]'' 



< OO. 



(5) 



(6) 



(7) 



Assumption 3. {Short range correlation) There exists > and ro > such 
that for A > and boxes G W'' oi 7/- {1 < k < n) with min^^i dist(y4fc, Ai) > 
r > tq and maxi<fc<„ dia.m{Ak) < r, 

< expi-rf^}, (8) 

where Ek{\) = {X^^JA,) < \}. 

Now we are ready to state our second result. 

Theorem 2. Suppose that AssumptionslJ\-lMhold and that there exists a reg- 
ularly varying function with index L > such that the integrated density of 
states associated with the operator in ([T]) admits the lower bound 

iV°(A) > exp{-0(l/A)(l + o(l))} as A ^ 0. (9) 

Then, there exists a constant Ci > 1 such that for any fixed x G W^, 

P-a.s. u^{t,x)>exY>{-Cit/tl){d\ogt){l + o{l))} as t ^ oo, (10) 

where ip is the asymptotic inverse of (p. 

Remark 1. The exponential behavior ([5]) and (Q of the integrated density of 
states is called the "Lifshitz tail effect" (cf. [12]) and is typical for the trapping 
Hamiltonian H^^. The index L is called "Lifshitz exponent". Using these ter- 
minologies, we can summarize our results as follows: if we have the Lifshitz tail 
effect with exponent L > 0, then \ogUi^{t,x) behaves like —t/{logt)^^^^°^^\ 



P( fl E,{X)j-nE,{X))F( fl E,{X)j 
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Finally we briefly comment on the relation to early studies on the quenched 
asymptotics of u^{t, x). We first give historical remarks. The first result in this 
direction has been obtained for the Brownian motion in the Poissonian traps 
by Sznitman [TS| (see also [IH]): 

(t,0) =exp{-ct/(logt)2/'^(l + o(l))} as t ^ oo, (11) 

with an explicit constant c > 0. The same asymptotics has also been proved 
for the discrete counterpart (the simple random walk in Bernoulli traps) by 
Antal [1]. These results are consistent to ours since in these cases, the Lifshitz 
exponent is known to be d/2 [131 [16]. Later, Biskup and Konig [2] considered 
the simple random walk in i.i.d. traps with more general distributions. A 
representative example in their framework is 

P(K;(0) <v)= exp as v^O (12) 

for some 7 G (0, 00). For such a model, they proved the quenched asymptotics 

P-a.s. M^(t,0) = exp{-xt/r(t)(l + 0(1))} as t 00 (13) 

with a constant x > described by a certain variational problem and a function 
r(t) = (logt)^/('^+^''')+°(^) (t — > 00) which is determined by a certain scaling 
assumption. It is remarkable that they also discussed the annealed asymptotics 
and as a consequence, the Lifshitz tail effect with the Lifshitz exponent {d + 
27) /2 was proved. Hence the relation we mentioned in Remark [T] has already 
appeared in this special class. 

Next, we comment on some technical points. The lower bound (Theorem 
[2]) is a slight modification of that of Theorem 4.5.1 in p. 196 of [19] and not 
genuinely new. We include it for the completeness and to use in an application 
given in Section On the other hand, the upper bound (Theorem [T]) contains 
some novelties. Besides the generality of the statement, our proof simplifies 
the existing arguments. To be more precise, our proof contains no localizing 
argument which all the proofs of above results rely on, see e.g. Lemma 4.6 in [2]. 
We will see in Section 14.11 that our result indeed gives a simple proof of the 
quenched asymptotics for the Brownian motion in the Poissonian obstacles. 

2 Proof of the upper bound 

We take k = 1/2 and a; = in the proof. The extension to general n and x are 
verbatim. Also, we give the proof only for the continuous setting. The proof of 
the discrete case follows by the same argument. We begin with the following 
general upper bounds for u^(t,x) in terms of the principal eigenvalue. 

Lemma 1. Under Assumption^ there exist constants 02,03 > such that 
uUt, 0) < C2(l + (A° ,{{-t, tryy/') exp { -A° ,{{-t, tr)t} + (14) 
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Proof. Let r denote the exit time of the process from {—t,tY. Then, by the 
reflection principle, we have 



(15) 



Now, (fT^ follows immediately from (3.1.9) in p. 93 of [19] under Assumption [H 

□ 

Due to this lemma, it suffices for ([6]) to obtain the almost sure lower bound 
for the principal eigenvalue ^^,i{{—t,tY). We use the following inequality for 
the integrated density of states 





exp < 




>-T>t 


+ Poir<t) 


< Eo 


exp < 


- j V^{X,)ds 


>;r>t 





A^°(A) > G N; \^,u{{-R,RY) < A}] 



> 



(16) 



which holds for any A > and R> 0. The first inequality is an easy applica- 
tion of the so-called "Dirichlet-Neumann bracketing" and can be found in [3], 
(VI. 15) in p. 311. Now, fix e > arbitrarily and let A = (1 — e)il){d\ogt)~^ and 
R = t. Then it follows from ffT6ll and (151) that 



(17) 



< (2t)'^exp {-m - e)-V(c^logt))(l + o(l))} 

^ 2d^d-d/{l~e)^il+o{l)) 



for some 6{e) > when t is sufficiently large. This right-hand side is summable 
along the sequence tk = and therefore Borel-Cantelli's lemma shows 



A°i((-4,4)') > {1 - e)^{dlogh)-' 



(18) 



except finitely many k, P-almost surely. We can extend this bound for all large 
t as follows: since ip{dlogt) is slowly varying in t, we have 



A°i((-t,t)'^)>A°i((-4,4)'^) 

> (l-e)^(dlogtfc)-i 

> {1 ~ 2e)ip{d\ogt)-^ 



(19) 



for tk-i < t < tk when k is sufficiently large. Combined with Lemma IH this 
proves the upper bound IQ. 
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3 Proof of the Lower bound 

We take k = 1/2 and x = again. Also, we only consider the continuous 
case. As in the proof of the upper bound, the principal eigenvalue plays a key 
role. Let us write X^{U) for the k-th smallest eigenvalue of — 1/2A in U with 
the Neumann boundary condition. Then we have another inequality for the 
integrated density of states 

N'^iX) < G N; \l,{{-R,Rr) < A}] 

- (2^^^^ ^ W,X^,{{-R,Rr) < X}F{Xl,{{-R,Rr) < X) (2°) 

<c,F{Xl,{{-R,Rr)<X), 

which holds for any A G (0, 1) and R> 0. The first inequality can be found in [3] 
again, (VI. 16) in p. 331, and the third one is a consequence of the classical Weyl 
asymptotics for the free Laplacian, see e.g. Proposition 2 in Section XIII. 15 
of [IS]. For arbitrary e > 0, let A = (1 + e)ijj{d\ogt)~^. Then, using ([20]) 
and (Q, we find 

F{Xl,{{-R,Ry)>X)<l-cl'N''{{l + e)^{d\ogt)-') 

< l-C4i(2t)-^/(i+^)^(i+°«) (21) 

< 1 - 

for some S{e) > when t is sufficiently large. 

Now we introduce some notations to proceed the proof. Let us fix a positive 
number 

M>i + | + i (22) 

and define 

J = {-t/{\ogtf',t/{\ogt)''Y n (logt)^^Z^ (23) 

a = i + {o,4j{diogty/^Y (zgj). (24) 

Note that miuj^j d(C*i, Cj) > diam(Cj) and both of them go to infinity as 
t oo. Therefore, by using fl2T]) and Assumption [3] recursively, we obtain 

P(A!!,i(Ci) > (1 + e)^{dlogty^ for all t G X) 

< n^W,i(^*) > (1 + e)V^(rflogt)-') +exp{-(logt)2} 

iez (25) 

< (1 _ t-d+5ie)y^iio,t)--^^' ^ exp{-(logt)2} 

< exp{-t'^(^)(logt)-'^^'} + exp{-(logt)2} 

for sufficiently large t. Since the right hand side is summable in t G N, Borel- 
Cantelli's lemma tells us that P- almost surely, 

there exists i G X such that X^^^id) < (1 + e)ip{d\ogty^ (26) 
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for all large t G M. The next lemma translates (l26l) to an upper bound for the 
Dirichlet eigenvalue: 

Lemma 2. There exists a constant ci > 1 such that F-almost surely, 

there exists i such that i{Ci) < ciip{d\ogt)~^ (27) 
for all large t. 

Proof. We choose d {i E I) for which X^^^iiCi) < (1 + e)ip{d\ogt)-\ This is 
possible for large t G N by f l26i) and then it also holds for all large t with slightly 
larger e by regularly varying property of ip. Let (pf denote the L^-normalized 
nonnegative eigenfunction corresponding to i(C'i) and d^Ci {i G X) the set 

{x G a;d(x,9C,) < eiP{d\ogtf'^}. (28) 

We further take a nonnegative function p G Cl{Ci) which satisfies 

p = 1 on \ d,Ci and || Vp||oo < 2e^^%l){d\ogtY^/^ . (29) 

Such a function can easily be constructed by a standard argument using mol- 
lifier. Substituting p0f to the variational formula for the principal eigenvalue, 
we obtain 

A°i(a) < / |V(p0f)P(x) + K.(x)(p0f)^(x)rfa;. (30) 

To bound the right hand side, we first use the uniform bound on eigenfunctions 
lloo < c^X^,i{.C^YI* (see e.g. (3.1.55) in p.l07 of ^) to see 

\\P4>1\\1> I ct>'^{xfdx>l-c,e. (31) 
Jc,\a,c^ 

Next, it is clear from (l29l) and the above uniform bound that 

/ |V(p0f)P(x) + VL(a:)(p0f)2(a:)da: 

<2 I \V<Y^\\x) + VYx)<p'^{xfdx + 2 I \VpWx)cp'^{xfdx ^^^^ 

< {2 + 9,CQt~^Y{d\ogt)~\ 

Taking e = (2cqY^ and plugging these bounds into fl30l) . the result follows. □ 

We also need the following almost sure upper bound. 
Lemma 3. Under Assumptions^ we have F-almost surely, 

sup V^x) < (3d log t)^/" (33) 

for sufficiently large t. 
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Proof. By Chebyshev's inequality, 



P( sup VUx) > (Srflogt)^/" ) 

< {AtfFf sup 14 (x) > {Sdlogty^A (34) 



sup exp{K;(x)°} 

.a;e[0,l]'' 



Since the last expression is summable in t G N, the claim follows by Borel- 
Cantelh's lemma and monotonicity of sup^g(_^^()d V^{x) in t. □ 

Now, we can finish the proof of the lower bound. We pick lv for which 
Lemma [2] and Lemma [3] holds. Then we can find a box Ci [i G X) satisfying 

X^JQ) < cMdhgt)-' (35) 

for sufficiently large t. Let (pf denote L^-normalized nonnegative eigenfunction 
associated with i(C'j). It is easy to see that there exists a box g+ [0, 1]*^ C Ct 
{q G Z"^) such that' 

Il0?lloo /" (Pf{x)dx> [ (Pfixfdxyl^idlogt)-". (36) 

We also know the following uniform upper bound: 

U?\\oo<c,X^,^iQf' (37) 

from (3.1.55) in [19]. Let us recall that the semigroup generated by has 
the kernel p^{s,x,y) under Assumption [T] (see Theorem B.7.1 in [Uj). We can 
bound this kernel from below by using the Dirichlet heat kernel P(-t,t)''is,x,y) 
in {—t,tY as follows: 

Puj{s, 0, y) > exp{ -s sup V^{x) U(_t,t)d(s, 0, y) 

xe(-t,ty > (38) 

>c-jS-'^l'^Q^^{-s(:id\ogtfl'^-c^\yfls) if |2/|<t/2, 

where the second inequality follows by Lemma [3] and a Gaussian lower bound 
for the Dirichlet heat kernel in [20j. Taking s = t/{\ogt)'^'^ and noting that 
|g| < 2s, we arrive at 

inf p^{s,0,y) > exp{-css/2} (39) 

j/Gg+[0,l]'' 

for sufficiently large t. 

Plugging fl35|) - fl39|) into an obvious inequality, we arrive at 

u^{t,0) = / puj{t,0,x)dx 

> / / Puj{s,0,y)pcu{t- s,y,x) ' dydx 

Jq+[0,1]'' \\<Pi Woo (40) 

>— j^exp{-A°i(C,)t-C8s/2} / (t>?ix)dx 
W'Pi Woo "'g+[0,l]'* 

> c#(logt)~'''/2exp{-cit/^(dlogt) - Css/2}, 
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where in the third hne, we have replaced p^, by the kernel of the semigroup gen- 
erated by with the Dirichlet boundary condition outside Cj. This completes 
the proof of the lower bound of Theorem [2] since s = t/(logt)*^ was chosen to 
be o{t /ijj {log t)). 

4 Examples 

We apply our results to two models in this section. The first is the Brownian 
motion in the Poissonian obstacles, where we see that our result recovers the 
correct upper bound. The second is the Brownian motion in a perturbed lattice 
traps introduced in [7], for which the quenched result is new. 

4.1 Poissonian obstacles 

Let us consider the standard Brownian motion [k = 1/2) killed by the random 
potential of the form 

V^{x) = J2w{x-uj,), (41) 

i 

where {u = J2i ^uji^^u) is a Poisson point process with intensity z/ > and W is 
a nonnegative, bounded, and compactly supported function. As is mentioned 
in Section [H Sznitman proved in [IB] the quenched asymptotics for this model: 

(t,0) =exp{-c(c/,z/)t/(logt)2/'^(l + o(l))} as t ^ oo, (42) 

where c{d, u) = Xd{i'uJd/dy^'^ with denoting the principal Dirichlet eigenvalue 
of -1/2A in B{0, 1) and cua = \B{0, 1)|. 

We can recover the upper bound by using classical Donsker-Varadhan's 
result [5j and Theorem [H Indeed, the above potential clearly satisfies Assump- 
tion [T] and the asymptotics of the integrated density of states 

logAr°(A) ~ -z/cUdA;J/'A-'^/2 as A^O (43) 

has been derived by Nakao [13] by applying an exponential Tauberian theorem 
to Donsker-Varadhan's asymptotics 

E[u^(t,0)] = exp|-g(d,z/)t^(l + o(l))} as t ^ oo (44) 
with ^ 

c{d, v) = -^{yuji) j . (45) 

Now an easy computation shows that the asymptotic inverse of the right hand 
side of dla]) is 

^{\) = X.^ucudr'/'X'/' (46) 
and then Theorem [1] proves the upper bound in (1421) . 

Remark 2. In this case, the lower bound given by Theorem [2]is not sharp as is 
obvious from the statement. (In the proof, we lose the precision in Lemma [2l) 
However, the lower bound can be complemented by a rather direct and simple 
argument in the Poissonian soft obstacles case, see [18]. So our argument 
simplifies the harder part. 
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4.2 Perturbed lattice traps 

In this subsection, we use our results to derive the quenched asymptotics for the 
model introduced in [7j. We consider the standard Brownian motion (k = 1/2) 
killed by the potential of the form 

V^{x) = Y,W{x-q-u,), (47) 



where ({t^glgez^, Pe) {0 > 0) is a collection of independent and identically 
distributed random vectors with density 

PeK G dx) = N{d, 9) exp{-|x|^}rfx (48) 

and is a nonnegative, bounded, and compactly supported function. The 
author has derived the annealed asymptotics for this model in [7] and also 
proved the following Lifshitz tail effect as a corollary: 



logiVD(A) 



-A-i-f (logi)-^ (rf = 2), 
(rf>3). 



where /(x) x^..^^. (y((x) means < Wmmix^^ f{x) / g{x) < limsup^^^ /(a;)/(7(a;) < 
oo. 

We can prove the quenched asymptotics from this result. 
Theorem 3. For any 6 > and x G Mf^, we have 

-t (logt)~2^(loglogt)"2^ ('^ = 2), 
logu^{t,x) Xt^oo { _ 2, (49) 

-t (logt) {d>3), 

with Fo-probability one. 

Proof. The Assumption [1] is clearly satisfied since K; is locally bounded almost 
surely. Hence the upper bound readily follows by computing the asymptotic 
inverse of (H9l) and using Theorem [H To use Theorem [2], we have to verify 
Assumptions [2] and [3l The former is rather easy and can be found in Lemma 11 
in [8]. The latter is verified as follows: we first fix ro > sufficiently large so 
that supp W C -6(0, ro/4). For r > ro and boxes {Ak}i<k<n as in Assumption[3l 
let us define events 

El = {for all qeZ'^ with d(g, ^i) < r/2, d(g + tOg, Ai) < 3r/4}, (50) 
E2 = {for all qeZ'^ with d(g, A^) > r/2, d{q + ujg, A^) > r/4}. (51) 

Then, 1(^1) and {A^ i(^fc)}2<A:<n are mutually independent on Ei fl E2 
thanks to our choice of rg. Therefore, the left hand side of ([8]) is bounded by 
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Fg{Ef)+Fe{E^). Let us denote the s-neighborhood of Ai by Ns{Ai). The first 
term is estimated as 

FeiEl) < Fe{\ujg\ > r/4 for some g G Z'^ n Nr/2{Ai)) 

<N{d,9)i^{qeZ'^nNr/2{Ai)} I exp{-\xf}dx (52) 

J\x\>r/4: 

< N{d,eyexp{-{r/8f} 



for large r, where we have used diam(Ai) < r in the last line. Next, we bound 
the second term Pe(£^2)- Using the distribution of oug, we have 

Fe{E!^) =Fe{q + u;ge Nr/iiAi) for some qeZ'^\ A^,/2(Ai)) 

<N{d,9) J2 [ exp{-\x - ql'^jdx 

,ezA7V./,(A0'^^'-/*(^^) (53) 

<N{d,ey Yl exp{-d{q,Nr/Mi)f}- 

q&d\N,/2{Ai) 

We can assume by shift invariance that Ai is centered at the origin. We divide 
the sum into two parts {|g| < r} and {|g| > r}. The former part of the sum is 
bounded by 

#{g G Z'^nfi(0,r)} sup exp{- d(g, iV,/4(A))^} 

qez^i\K/2iAi) (54) 

< Cior'^exp{-(r/4)^}. 

For the latter part, we use the fact that Nr/4{Ai) C i?(0,3r/4), which follows 
from the assumption diam(74i) < r. By using this fact, we find 

d(g, Nr/Mi)) > \q\ - ^r/A > |g|/4 for \q\ > r (55) 

and therefore 

exp{-d(g,iV,/4(A))^}< exp{-|g/4r}. (56) 

q&<i\N^I2{Ai),\q\>r q&'i,\q\>r 

It is not difficult to see that this right hand side is bounded by exp{ — (r/8)^} 
for sufficiently large r. Combining all the estimates, we arrive at 

Pe(^D + Fg{E',) < N{d, ey{2 + cio/) exp{-(r/8)^} (57) 

for large r, which verifies Assumption [3l □ 
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